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ABSTRACT 


By separating the microscopic equations of nation for 
the particle density and the fields in a plasma into 
ensemble average and fluctuating parts we have been able 
to obtain an unambiguous description of scattering of an 
electromagnetic > wave from density fluctuations. The 
method described can be generalized to include all linear 
scattering processes in a plasma involving the inter- 
action of transverse and longitudinal waves with fluccu- 


at ions 



DRAFT 


I. INTRODUCTION 


In an earlier paper , the author had attempted an un- 
ambiguous formulation of the cross section describing 
fluctuation scattering and absorption in a space-dispersive 
collisionless plasma. In the process, an ansatz was in- 
voked concerning the introduction of fluctuations in the 
Boltzmann- Vlas sov (B-V) equation for a plasma. As is well 
known, the B-V equations are ensemble average equations for 
the particles and fields and as a consequence do not re- 
flect the existence of fluctuations in a plasma; hence, the 
need for the aforementioned ansatz. In the present paper 
the same problem will be discussed in a more rigorous manner 
by introducing the fluctuations in the microscopic equations 
of motion for the particles in the plasma. 

The need to introduce statistical mechanical methods 
for the systematic derivation of the macroscopic Maxwell's 
equations in a medium and of fluctuation scattering has been 
recognized many years ago(^#3). It is however, relatively 
recently that Fixman^ 4 ^ attempted the description of the 
problem by simultaneous use of microscopic and macroscopic 
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equations for the electromagnetic field. The implicit 
assumptions in Fixman's work were that the microscopic 
equations truly represent the complete picture of the 
microscopic events in the interaction of fields with 
granulated matter and that the macroscopic equations 
truly represent the ensemble average properties of the 
same interaction. Simultaneous manipulation of the 
microscopic and ensemble average fields has yielded a 
sensible, if not completly* satisfactory^ approach to 
the problem. More recently Mazur (5, 6) 'has also dis- 
cussed the statistical mechanics of the electromagnetic 
properties of matter and has given a molecular theory 
of light scattering that closely follows Fixman's steps. 

In what follows, we will develop the theory of fluctu- 
ation scattering in a collisionless plasma in a more 
natural way starting from the microscopic equations for 
particles and fields. It is believed that this theory 
can adequately and without ambiguity describe all linear 
scattering processes in a plasma (transverse waves from 
fluctuations, longitudinal waves from fluctuations, fluctu- 
ations from fluctuations) whose particle correlations higher 
than the second can be neglected. 
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II. THE MICROSCOPIC EQUATIONS FOR A PLASMA 

For the microscopic description of the plasma, 
we follow the method initiated by Klimontovich ) . He 
introduces the random function of the microscopic 
particle density for the species in a multicom- 

ponent plasma as 

Nfifih = (1) 

j J 

to define the number of particles in phase space apaq 
around p and q at time t.p and q are canonical variables, 
Wj, is the number of particles of the •'r-th species. From 
the Hamiltonian of the system (particles and fields) 
Klimontovich then derives the following exact equations 
of motion for (f^i) and of the electromagnetic field 
c/ t where A(tj'i) is the vector potential and 
the scalar potential 
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where 


* * l) "Z e - ( rfTT ld ^' 


y ~ If i'i (4) 

As the discussion will not include an external magnetic 
field the term — (j*x has been omitted in Eq. 2. 

It should be noted that there 'arc as many equations 2 
as there are species in the system under consideration. 

Our objective is to begin with the foregdj^ng equa- 
tions and derive macroscopic (ensemble average) equa- 
tions for the electromagnetic field that will describe 

the index of refraction and fluctuation scattering 

( 1 ) 

processes. We make the following assumptions' . 

(a) The microscopic quantity can be 

written as a sum of its ensemble average mlue 

and a fluctuating part 

Wj'iyiji) -r having a zero ensemble average: 

< (?&/>- < /(pf*'-' 

(5) 


(b) The microscopic electromagnt ; c fie] a 
can be similarly separated in ensemble average n a 
fluctuating parts 

Amu - A 

- f (<:M . V' 
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with 

<A6jt;>- A'*’(<iU ^ <A'\u> -O 
<4>MU> -* ^ - 0 . 

The separation of the fields in ensemble average and 
fluctuating parts is significant. The ensemble average 
parts describe the average properties in a plasma such 
as propagation, index of refraction, etc., i.e., the 
macroscopic properties. In this description the plasma 
has lost the granularity implied by Eq.-l; it is a con- 
tinuum. The fluctuating parts describe the tering 

due to the incomplete extinction of the microscopic prop- 
agating fields and, on the average, are zero, inis 
statement corresponds to Yvon's prescription for the 
evaluation of fluctuation scattering. It should be noted, 
that when only a transverse propagating ensemble average 
wave exists^ ) itself is zero. 

We now introduce (5) and (6) in equations (2) and 
(3) . For equation (3) we then have 


( T - ?£)a‘v> - fI^) + \m m w 


H I 

1 
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On taking ensemble averages in Eq. (7) only terms in 

survive; subtracting the resulting en- 
semble average equation from Eq. 7 we also obtain an 
equation for the fluctuating quantities A , A/^ 

We then have the system of equations 







jrj . 

4> qO 




(9) 


where <T«0,J . We now turn our attention to the equation 
of motion for . We introduce the operator 


o lipj+J 


i t - \7~ 


and henceforth write 
Eq. 2 then takes the form 


X(p +- e r 

-+ 

i £(p jt) ^ t . 1 '(ijt '** V- /-y# Vi L 


( 10 ) 



- 7 - 


The first step "here is to take the ensemble average of 
Eq. 10. Clearly, what remains of Eq. (10) is 

Because of Eq. 8, of Eq. 11 represents quantities 

of the form < /¥"($'*) > - M 

i.e., two-point, one-time correlations of the fluctuations. 
In the derivation of the celebrated Vlassov equation these 
terms are neglected thus breaking the coupling with the 
hierarchy of equations for the higher-order correlations. 
Parenthetically we note here that the solution of the 
Vlassov equation (Eq. 11 with ^•••✓ > ~0) is obtained by 
writing t'l^C jp<j c ) - - f (ptjLj and by neglecting 

the non-linear terms V ^'(Lt j • These simplifications 

lead to the calculation of the space-time dispersive 
dielectric constants 

We now subtract Eq. 11 from Eq. 9 and obtain 
tCy £ t0 (<(' t)- V. -i- - 'U / • *y f 


* -f«j , 

+ S t " 




(12) 


where the double b^r implies that for a quantity F , F a? F -O r >‘» 



- 8 - 


( 

In the absence of ensemble average fields in the system j 

and on the assumption that two-point one-time correlations | 

can be neglected (i.e. r-0 ) eq. 12 reduces to j 

t l \t) =--0 (12a) | 

r ^ ^ i 

where £ 0,, c ~ J- Cl!'* and where j>/> Q~ are the * 

non-driven (naturally occuring) field fluctuations in the 

plasma. This equation will be useful later in the evaluation 

» f «) 

of correlations of fluctuations in 

Equations 8, 9, 11, and 12 are exact equations con- 
taining both the ensemble average and the fluctuations of 
particles and fields. They provide the essential quantities 
..ceded to describe propagation of transverse and longitudinal 
waves, scattering of fluctuations from fluctuations, etc. 

For the purposes of the present paper we shall make the 
following simplifications: 

(a) The ensemble average propagating transverse 1 

electromagnetic field interacts only little with the heavy 

ions in the system. It is, therefore, sensible to keep 
one term in the summation of the right hand side of Eq. 8 
representing the electron current. 

(b) To obtain results amenable to calculation we 

, . i/(®> 

will use the standard perturbation expansion «■ 
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with ^ and nonlinear terms such as W‘/« 


P i> M 

will be neglected. 


III. THE HIERARCHY OF EQUATIONS FOR THE FIELDS 


It becomes expedient at this point to transform the 
equations derived earlier from to 4ft0-space . With 

r - ■ L. . icjt' 

= \k^ 1 fA% 


we have for Eq. 8 (neglecting all particles except 
electrons) 




•Hi 


' - r j 


(13) 


For *Eq. 9 




y 


*We have = V i>(j) • **• 

Clearly then after a Fourier tr ansjformation 

- (4 , ■= - 4Wx~ ^ . However , 1 Nf&djt 

hence, the result of Eq. 14. 
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For Eq. 11 (with *>VN®(p) and /(p*«0 — tU + iip/'Wf ) : 

f i 1 # r 

{ pis vo) Ai^ (p-foo) + ? V&o) • V- 1^/ ( p ) 


4 e. 


< £ 


and for Eq. 12 (again 




(15) 


): 


/ip-U w' w (pU + E*««o>. v, fJ> ) 


r 

4 €. 


4 €. 


fr'iHnO'? w>-w') •liV 

p r* ^ ^ £?r)* 


0 

(16) 


We can now solve for U}^ and ^ and introduce the 

HoJ/ (■ 

results in Eq. 13. The contributions of t Vp-/^ iPJ 

from Eq. 15 and of e (,> $«».v. wJ,°V ) from Eq. 16 convert 
the free space operators of the left hand side of Eq. 13 
to operators acting on the ensemble average fields in a 
space-dispersive medium. Explicit ely, and introducing the 
transverse and longitudinal dielectric constants £.(-fcu)) and 
, we obtain 


-II- 




X,(L) A“\U) + = 


_ 4£ 
c 


w MM 

and similarly 

*X u (ii») ^(L) 

_ ;'4irjif f «*V dp£ 


-•deSS* 

(17) 


I' 


, f£ w (lv)-vr(p, U>-«o 

£{pfa) L p P 


where 




(18) 


$ r (&wo) =-4- G^c)*" Gj&ti) j XjArt - & (obfc ) 6j ^uJ) 


r *r 4rre] p 

6 r (^u)) =* ( + .A 


^ 0)- M * * f~ 




I 


1 4*& 

- 4 * a* 


,(oa x d tr 


l\JJJ% v y r , .. 

^ t*)* ■■?■'• sL/‘* -t xi ^ 


Equations 17 and 18 form the desired set of coupled equations 
for the ensemble average (A^4' / and the fluctuating 
J fields. Except for the linearization, these 
equations are exact equations and by themselves form a 
hierarchy of equations for higher order correlations of 
particles and fields." All correlations above the second 
will henceforth be neglected by omitting the F term in Eq. 18. 

To fix ideas on how to proceed from this point we will 
assume that the only ensemble average field in the system 
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is a transverse wave such as would be excited by an in- 
cident electromagnetic wave. Then so that 


1 ‘ 


n 


- 1 


C /' Hf t 


tl),l 






(19) 





f\jr- 


gj %L mtwwcJ-iM 


( 20 ) 

.0 ff 

Multiplying both terms of Eq. 20 by 1/t and and using 

the fact that J>/\ or -k'A ~0 we can separate the transverse 
from the longitudinal components. We then find 



j'4r\e^ 
v c Jm t 




f r 


')4>_ 

0«f? .u>, 

(21) 


» 

W h!AA 


whence 


r, (p,M 


fc). 






i 



( 22 ) 


( 23 ) 
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l(0> 70; .(,) 

with expressions now available for f\ y n and q> we 

can now proceed to evaluate the scattering cross sections. 


IV. THE SCATTERING CROSS SECTION 

\ 


As was indicated in Ref. 1 a possible method for the 
evaluation of scattering cross sections consists of in- 
troducing the values of A**/ $ ^ from Eq. 21 into Eq. 20 

and then evaluating the changes this brings on the imagin- 
ary part of This would give a direct and unambiguous 

measure of the cross sections. A difficulty, however, oc- 
curs here because of the presence of terms having the form 

, necessitating a rather 

intricate calculation. Although these terms have been 
evaluated in Ref. 7 we shall presently resort to some 
simplifications in the momentum space integratic which 
will hopefully make the results less opaque. 

Take for instance the quantity appearing in Eq. 21: 

i,« UK- (ip k.U'A--* 

j£.(pjkjA) r ' 

* 

Using successive vector identities we find that 


r.. ((&■})* 




s 'flu shorthand 'notation usui itctw is . 
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Converting the integral over the divergence to a line 
integral and noting that for p-y«O / ^-»0 we readily obtain 

T. _ p p w;v4 f jA ( “'.^^j 

Assuming further that u 2 > 4 .p/n f i-e. that £~>-ild we obtain 

T . - i U N%, KV ■» 

0 »>£-p/>V J r u) J* " 

- i ) t« w w < m; 


66 * 


T <0 


where vj is the current fluctuation. 

Note that the continuity equation would imply 

(K-.tj - K’Jp (Kj£l )-~ 0 thus X^O • In a similar fashion 
we find that 

Integrations in p of Eq. 19 can also be simplified. We 


find 




V'/ *1 


-A(fM) 




v - 


r 

fe-p/rtl 


- 1 A"’(-;y; .£'$a) -v £ v.4 w (W >OT"fca> 



f 


A f ' ) 
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in Eq. 19 


It now becomes obvious that inserting t f j\ 
leads to terms $ 

etc. Restricting ourselves only to density-density 
correlations we then find for Eq. 19 when &->> i'-p/ui 

W AW - lt)^1 IS Si 


VJ fcJir/ 
o' , a" 


pi 

[woi u 


J" 1 < 

^ c 




(24) 


However (1) 




'^'* v "v’ v **/’y 

so that finally 


A 


(26) 


where 




C*2JTJV 




> (27) 


pttu; 


1- (a.-t*)4 


- + 


a r» > 

c 4 (-».* , 


(28) 


-l*k) e r (&4) 


~ n cd u -&' 4 c. (£W> 
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A ,(o) * «/#) 

and 2, is a unit rector for the rector potential: A 
The scattering cross section , can then be evaluated *from 
■X'vyj •[ (■&&)} • The fi rst term in Eq. 28 describes the 
scattering of a transverse wave on density fluctuations to 
produce a transverse wave; the second term describes the 
scattering of a transverse wave on fluctuations to produce 
a longitudinal wave. If all the terms in Eqs. 17 and 18 
(except F ) were kept, additional linear scattering processes 
could be described. 

the cross section GY-few)), gives a generalized r ersion 
of the Prins-Zernike formula describing scattering from 
the covariance of density fluctuations, for it also includes 
"inelastic" scattering processes . 


T(4<o) — 2 A/ 0 ~' (ui/c) V/ i 

where Oj a (*«))=-£ r (k<j) and where N, is the particle density. 

The evaluation of T and the consequences of the second term 
of F/Mdfj involving ) are discussed at some length <t>n 


Ref. 1. 
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¥ 


CONCLUSIONS 


By separating the microscopic equations for particles 
and fields in a plasma into ensemble average and fluctu- 
ating parts we have been able to derive a generalized 
version of the Pr ins-Zernike formula for scattering of an 
electromagnetic wave from the covariance of density 
fluctuations in a plasma. The procedure leading to the 
desired results is shown schematically in Table I. The 
method developed can be extended to cover all linear 
scattering processes involving the interaction of trans- 
verse and longitudinal waves with fluctuations in a plasma. 
This method does not seem to be restricted to the des- 
cription of small fluctuations, for nowhere did we have 
to resort to "small signal" analysis. 


*• ► 
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